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ABSTRACT 

t 

/ 

J 

A  simultaneous  numerical  and  experimental  investigation  is  undertaken  to  assess  three  dimen¬ 
sional  effects  and  HRR  dominance  near  a  crack  front  in  a  ductile  3-point  bend  specimen.  In 

/ 

parallel  to  the  3-D  numerical  calculations,  a  plane^strain  and  a  plane-stress  analysis  of  the  same 
in-plane  specimen  geometry  is  performed  to  obtain  upper  and  lower  bounds  for  the  3-D  calcula¬ 
tion.  The  radial,  angular  and  thickness  variation  of  the  stresses  and  displacements  are  studied  in 
great  detail  from  contained  yielding,  to  fully  plastic  conditions.  The  results  indicate  that  the 
plane  strain  HRR  field  prevails  in  the  interior  of  the  specimen  very  near  the  crack  front  even  for 
moderate  extents  of  yielding.  On  the  other  hand,  for  distances  from  the  crack  tip  exceeding 
about  half  a  specimen  thickness,  plane  stress  conditions  are  approached. 

The  calculations  presented  here  model  a  series  of  laboratory  experiments  involving  three 
independent  experimental  techniques.  Details  regarding  the  experiments  and  comparisons  of  the 
experimental  measurements  with  numercial  calculations  and  theory  are  emphasized  in  Part  II  of 
this  work.  >1  u  >  t  • 


1.  INTRODUCTION 

Fracture  in  structural  components  of  practical  dimensions  which  are  made  of  low-to- 
intermediate  strength  metals  is  often  accompanied  by  extensive  plastic  deformation  prior  to 
crack  initiation.  Hence,  one  is  compelled  to  seek  elastic-plastic  solutions  to  crack  problems  to 
treat  such  situations.  The  HRR  solution,  (Hutchinson,  1968,  and  Rice,  1968)  describes  the 
asymptotic  stress  and  strain  variations  near  a  monotomically  loaded  stationary  crack  tip  in  a 
material  characterized  by  the  deformation  plasticity  theory  under  two  dimensional  (plane  strain, 
or  plane  stress  or  anti -plane  shear)  conditions.  The  J  integral  (Rice,  1968)  which  may  be  inter¬ 
preted  as  the  energy  release  rate  for  a  material  obeying  the  deformation  plasticity  theory  also 
provides  a  measure  of  the  intensity  of  the  near-tip  fields  in  the  HRR  solution. 

As  emphasized  by  Hutchinson  (1983),  it  is  imperative  that  the  HRR  solution  should  dom¬ 
inate  over  a  size  scale  which  is  large  as  compared  with  the  region  near  the  tip  wherein  finite 
strain  effects  are  important.  The  above  requirement  ensures  that  J  can  be  safely  used  to 
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characterize  the  initiation  of  crack  growth  in  ductile  materials. 

It  has  recently  been  demonstrated,  on  the  basis  of  plane-strain  calculations,  that  finite  strain 
effects  become  important  only  for  distances  from  the  tip  of  the  order  of  2  or  3  times  the  crack 
tip  opening  displacement. 

The  issue  of  /-dominance  has  been  examined  by  many  investigators  in  a  variety  of  situa¬ 
tions,  ranging  from  contained  yielding  to  fully  plastic  conditions  under  Mode  I  plane  strain 
(Shih  and  German,  1981  and  McMeeking  and  Parks,  1979).  As  an  outcome  of  these  studies,  it 
has  become  possible  to  specify  proper  dimensions  for  some  commonly  used  specimens  in  frac¬ 
ture  toughness  testing,  so  that  /-dominance  could  be  ensured.  As  noted  by  Hutchinson  (1983), 
analogous  specifications  for  plane  stress  have  not  as  yet  been  obtained.  Nevertheless,  the  issue 
of  HRR  dominance  has  partly  been  addressed  by  Narasimhan  and  Rosakis  (1988),  by  means  of 
a  detailed,  numerical  calculation  which  modeled  Mode  I,  plane  stress  cracks  under  small  scale 
yield  conditions. 

All  the  analyses  mentioned  in  the  foregoing  section  are  based  on  the  assumption  of  two 
dimensional  plane  strain  or  plane  stress  conditions.  No  analytical  asymptotic  solution  exists 
near  three  dimensional  crack  fronts  in  elastic-plastic  materials.  In  view  of  the  above,  it  is 
important  to  examine  the  regions  near  the  crack  front  in  which  the  plane  stress  and  plane  strain 
HRR  solutions  provide  a  good  approximation  to  the  actual  three-dimensional  fields.  It  is  gen¬ 
erally  assumed  (see  for  example  Li,  Shih,  and  Needleman,  1985)  that  plane  strain  conditions 
prevail  as  the  crack  front  is  approached,  at  least  in  the  interior  of  the  specimen.  But  it  appears 
that  no  detailed  investigation  of  the  region  of  dominance  of  any  singular  elastic-plastic  field  has 
been  carried  out  under  three-dimensional  conditions.  This  issue  is  pivotal  to  the  definition  of  a 
local  energy  release  rate  (pointwise  along  the  crack  front),  see  for  example  (Li,  Shih  and  Nee¬ 
dleman,  1985,  and  Nakamura,  Shih  and  Freund,  1986a,  and  1986b).  and  its  interpretation  as  an 
amplitude  factor  of  the  near-tip  fields. 

Surface  integral  generalizations  of  the  line  integral  expression  of  /  to  three-dimensional 
crack  problems  are  available  in  the  literature  (see  for  example  Budiansky  and  Rice,  1973). 
Nevertheless,  the  definition  of  a  local  (pointwise)  energy  release  rate  and  its  numerical  compu¬ 
tation  has  been  the  subject  of  much  investigation  (Li,  Shih  and  Needleman,  1985,  Shih,  Moran 
and  Nakamura,  1986,  Nakamura,  Shih  and  Freund,  1986a,  Delorenzi  and  Shih,  1983,  Wellman, 
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Rolfe,  and  Dodds,  1987,  and  Dodds,  Carpenter  and  Sorem,  1987). 

In  Delorenzi  and  Shih,  (1983),  Wellman,  Rolfe  and  Dodds  (1987),  Dodds,  Carpenter  and 
Sorem  (1987),  and  Nicholas  and  Luxmoore  (1987),  fracture  parameters  such  as  J  and  crack 
opening  displacement  are  computed  using  three  dimensional  numerical  solutions  for  compact 
tension  and  three  point  bend  specimens  and  are  compared  with  experimental  results.  The  varia¬ 
tion  of  three  fracture  parameters  and  of  the  stresses  through  the  thickness  are  also  examined  in 
Delorenzi  and  Shih,  (1983)  Wellman,  Rolfe  and  Dodds  (1987),  Dodds,  Carpenter  and  Sorem, 
(1987),  and  Nicholas  and  Luxmoore  (1987). 

The  purpose  of  the  present  investigation  is  twofold.  Firstly,  numerical  and  experimental 
methods  are  simultaneously  used  to  assess  the  effects  of  three  dimensionality  near  a  crack  front 
in  a  ductile  three  point  bend  specimen.  The  experiments  are  based  on  the  optical  methods  of 
caustics  and  interferometry  and  are  described  in  detail  in  Part  II  of  this  investigation.  (Zehnder 
and  Rosakis,  1988).  Some  results  pertaining  to  the  above  issues  and  also  to  the  variation  of  J 
and  crack  opening  displacement  will  be  given  in  Sec.  4. 

The  second  issue  that  is  addressed  here  (see  Sec.  5),  is  the  extent  of  regions  over  which 
plane  strain  and  plane  stress  analyses  provide  a  good  approximation  to  the  actual  three  dimen¬ 
sional  fields.  To  this  end,  detailed  three-dimensional  and  two-dimensional  numerical  calcula¬ 
tions  have  been  performed  (see  Sec.  2).  These  calculations  model  the  monotonic  loading  of  a 
ductile  three  point  bend  specimen  from  contained  yielding  to  fully  plastic  conditions.  The 
radial  and  angular  variation  of  the  stresses  from  these  analyses  will  be  compared  with  the  HRR 
solution  in  the  spirit  of  Shih  and  German,  (1981),  McMeeking  and  Parks  (1979)  and 
Narasimhan  and  Rosakis  (1988),  with  the  view  of  assessing  J-dominance.  Some  preliminary 
results  pertaining  to  the  above  issues  have  been  discussed  by  Narasimhan,  Rosakis  and  Zehnder, 
(1988). 

2.  NUMERICAL  ANALYSIS 

The  three-point  bend  specimen  shown  in  Fig.  la,  was  modelled  using  three  dimensional 
finite  elements.  The  geometry  of  the  mesh  is  shown  in  Fig.  lb.  Due  to  symmetry,  only  one 
quarter  of  the  specimen  was  simulated,  with  appropriate  boundary  conditions  imposed  on  the 
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planes  of  symmetry.  The  mesh  consists  of  5  layers  of  elements  (through  half  the  thickness  of 
the  plate).  Each  layer  is  composed  of  420  eight-noded  brick  elements.  The  layer  interfaces  are 
located  at  x$/h  =  0,  0.15, 0.275, 0.375,  0.45  and  0.5.  It  should  be  noted  that  the  layers  become 
thinner  as  the  free  surface  is  approached  in  order  to  model  the  comer  singularity.  The  layer  of 
elements  adjoining  the  mid-surface  of  the  plate  is  three  limes  as  thick  as  the  one  near  the  free 
surface. 

The  crack  tip  elements  are  formed  by  collapsing  the  eight  noded  brick  elements  to  triangular 
prisms.  The  details  of  the  in-plane  mesh  near  the  crack  tip  is  shown  in  Fig.  lc.  This  fine  mesh 
is  composed  of  a  circular  region  whose  radius  is  equal  to  the  thickness  of  the  plate.  There  are 
ten  circular  rings  in  this  region,  which  are  divided  by  19  rays  at  equal  angular  intervals  of  10°. 
The  radius  of  the  innermost  ring  of  elements  (nearest  to  the  tip)  is  1/100  of  the  plate  thickness. 
This  detailed  mesh  is  expected  to  provide  sufficient  resolution  near  the  crack  tip,  so  that  a  care¬ 
ful  study  of  the  asymptotic  stress  and  strain  variations  may  be  performed. 

Two  dimensional  plane  strain  and  plane  stress  calculations  were  also  performed  with  the 
same  in-plane  mesh  (Fig.  lb)  to  obtain  upper  and  lower  bounds  for  the  3  dimensional  calcula¬ 
tions. 


A  small  strain,  J2  incremental  plasticity  theory  was  used.  The  material  was  assumed  to  be 
homogeneous,  isotropic  and  to  obey  the  Huber-Von  Mises  yield  criterion.  The  response  of  the 
material  in  uniaxial  tension  was  characterized  by  a  piecewise  power  hardening  law  of  the  form: 


—  =  J  C '  °°  a-°o 
to  |  {o/aey  a>o0 


(2.1) 


with  a  hardening  exponent  n  =  22  and  yield  stress  o0  =  1030  Mpa.  These  values  were  chosen 
to  match  the  constitutive  properties  of  the  particular  heat  treatment  of  4340  carbon  steel  used  in 
the  experiments.  These  experiments  are  described  in  detail  in  part  II  of  this  investigatioa 

The  B-bar  method  proposed  by  Hughes  (1980)  and  later  modified  by  Nakamura,  Shih  and 
Freund  (1986)  was  used  to  relieve  artificial  mesh  locking  effects  that  occur  under  nearly 
incompressible  conditions  (in  the  fully  plastic  range)  under  plane  strain  and  in  three  dimensions. 
Nakamura,  Shih  and  Freund  (1986a)  and  (1986b)  suggested  that  the  B-bar  method  should  be 
used  with  a  stabilization  procedure  akin  to  Belytschko  and  Tsay  (1983)  to  safeguard  against  a 
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spurious  pressure  mode.  As  recommended  in  Nakamura,  Shih  and  Freund  (1986),  a  stabiliza¬ 
tion  parameter  of  0.05  was  used  in  the  present  calculations. 

The  stress  computations  were  performed  using  the  tangential  predictor-radial  return  method 
with  subincrementation  (see  appendix  of  Narasimhan  and  Rosakis  1988).  The  nonlinear  finite 
element  equilibrium  equations  were  solved  incrementally  using  a  quasi-Newton  method  (BFGS) 
(Matthies  and  Strang,  1979).  For  the  three-dimensional  case,  this  method  provided  substantial 
savings  as  compared  with  the  full  Newton  method.  The  convergence  to  an  equilibrium  solution 
was  quite  rapid  in  the  early  part  of  the  load  history,  requiring  typically  3  or  4  iterations  in  each 
load  step.  But  it  was  quite  slow  as  fully  plastic  conditions  (given  by  the  rigid  plastic  limit  load) 
were  approached.  As  suggested  in  Matthies  and  Strang  (1979),  a  BFGS  update  was  not  per¬ 
formed  whenever  the  condition  number  of  the  updating  matrix  became  very  large  (greater  than 
10 5).  This  happened  quite  frequently  when  the  limit  load  was  approached.  Also,  smaller  load 
increments  were  used  near  the  limit  load. 


3.  DESCRIPTION  OF  EXPERIMENTS 

As  described  in  detail  in  part  II,  the  experiments  were  performed  on  three  point  bend  speci¬ 
mens  of  a  ductile  4340  steel  with  dimensions  given  in  Figure  la.  Two  specimens  were  tested. 
In  both  cases,  the  applied  load  and  the  load  point  displacement  were  measured  and  recorded 
during  the  experiments.  In  order  to  provide  additional  comparison  between  the  numerical  cal¬ 
culation  and  the  experiments,  a  measurement  of  the  out  of  plane  displacement  fields  of  the 
specimen  surface  was  performed  using  Twyman-Green  interferometry. 


4.  RESULTS 

4.1  Load-Displacement  Curves 

The  variation  of  load  versus  load  point  displacement  is  shown  in  Fig.  2  from  the  three 
dimensional  and  two  dimensional  numerical  calculations.  For  comparison,  the  experimental 
load-displacement  records  for  two  specimens  are  shown  in  the  figure  up  to  the  point  of  fracture 
initiation.  The  plane  stress  and  plane  strain  limit  loads  (for  a  rigid-perfectly  plastic  material) 
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are  also  indicated  in  the  figure.  The  curve  obtained  from  the  3-D  analysis  lies  in-between  the 
curves  from  plane  stress  and  plane  strain  calculations,  and  is  much  closer  to  plane  stress.  Also, 
very  good  agreement  can  be  seen  between  the  experiments  and  the  3-D  computation. 

4.2  The  J  Integral 

The  J  integral  defined  by  Rice  (1968)  for  a  two  dimensional  cracked  body  is  path  indepen¬ 
dent  for  a  material  obeying  the  small  strain  deformation  theory  of  plasticity.  Under  the  above 
conditions,  this  integral  may  be  inteipreted  as  the  energy  release  rate  and  is  also  a  characteriz¬ 
ing  parameter  of  the  asymptotic  crack  tip  fields  (Hutchinson,  1968  and  Rice  and  Rosengren, 
1968). 

For  a  three-dimensional  crack,  the  J  integral  is  defined  over  a  cylindrical  surface  (Budian- 
sky  and  Rice  1973,  and  Broberg  1987)  surrounding  the  entire  crack  front.  The  numerical  com¬ 
putation  of  the  above  surface  integral  in  a  3-D  finite  element  analysis  is  very  difficult.  How¬ 
ever,  the  surface  integral  may  also  be  interpreted  as  the  energy  released  due  to  a  unit  (normal) 
virtual  crack  extension  along  the  entire  crack  front.  An  average  value,  denoted  here  by  J ,  can 
then  be  obtained  by  dividing  the  above  quantity  by  the  length  of  the  crack  front.  (Nakamura, 
Shih  and  Freund,  1986  and  Delorenzi  and  Shih,  1983). 

If  5/(j)  denotes  the  crack  advance  at  a  point  s  along  a  three  dimensional  crack  front  (in  a 
direction  normal  to  it)  and  d s  the  elemental  arc  length  along  the  crack  front,  a  local  quantity  J 
Qc3  /  h)  has  been  defined  by  Li,  Shih  and  Needleman  (1985)  and  Shih,  Moran  and  Nakamura, 
(1986): 

~  J  /  (jc3//i)  5/(jc3  /  h)dxi  =  -  5 n 

"  -44 

Here  811  is  the  change  in  potential  energy  of  the  body.  It  has  been  (heuristically)  argued,  see 
for  example  Li,  Shih  and  Needleman  (1985)  and  Shih,  Moran  and  Nakamura  (1986)  that  j(s) 
plays  a  role  of  characterizing  parameter  (or  amplitude  factor)  of  the  crack  tip  fields  provided 
plane  strain  conditions  prevail  through  the  thickness  as  the  crack  front  is  approached  and  assum¬ 
ing  that  proportional  loading  is  pointwise  enforced.  As  noted  earlier,  this  argument  can  be 
accepted  only  after  a  complete  investigation  of  the  three-dimensional  near-tip  fields  validates 
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the  above  assumptions. 

A  domain  integral  representation  of  J  and  J(s)  has  been  discussed  by  Li,  Shih  and  Needle- 
man  (1985),  Shih,  Moran  and  Nakamura,  (1986),  Nakamura,  Shih  and  Freund  (1986a),  and 
Nakamura,  Shih  and  Freund  (1986b).  This  method  can  be  easily  implemented  in  a  finite  ele¬ 
ment  program  and  has  also  been  found  to  give  very  accurate  results.  In  particular,  it  is  very 
effective  in  computing  the  near-tip  value  of  J.  This  method  was  employed  in  the  present 
analysis  to  calculate  J  from  the  numerical  solution. 

The  variation  of  J  versus  the  applied  load  is  shown  in  Fig.  3a.  Also  shown  in  the  figure  is 
the  J  versus  applied  load  record  calculated  from  the  experimental  result  on  the  basis  of  meas¬ 
ured  load-displacement  curves  (see  Part  II).  In  Fig.  3b,  J  is  plotted  against  the  load  point  dis¬ 
placement  It  should  be  observed  that  the  curve  corresponding  to  the  3-D  solution  is  again 
bounded  by  the  plane  strain  and  plane  stress  curves.  For  small  load  levels,  the  three  analyses 
show  results  which  are  not  very  different  from  each  other.  These  values  also  agree  very  well 
with  J  calculated  on  the  basis  of  small  scale  yielding.  Also,  J  from  the  3-D  analysis  is  close  to 
the  plane  stress  solution  for  the  entire  load  range,  except  very  near  the  plane  stress  limit  load. 

It  is  worth  emphasizing  here  that  the  global  3-D  specimen  behavior  such  as  the  load  dis¬ 
placement  and  the  average  J  integral  records  (see  Fig.  2,  3)  are  closer  to  the  plane  stress  rather 
than  plane  strain  approximations.  This  is  not  surprising  given  that  the  ratio  of  thickness  to  the 
in-plane  dimensions  of  the  specimen  were  relatively  small. 

Nevertheless  as  will  be  shown  in  the  following  section,  at  the  vicinity  of  the  crack  tip,  there 
is  a  strong  variation  of  fields  through  the  thickness  which  indicates  the  need  for  a  3-D  analysis 
for  the  study  of  crack  initiation.  In  Fig.  4,  J  normalized  by  J  is  plotted  against  normalized 
distance,  x3  /  h ,  along  the  crack  front  for  different  ratios  of  the  applied  load  P  to  the  plane 
stress  limit  load  P0.  The  variation  for  the  elastic  case  is  also  shown  in  the  figure.  It  is  seen 
that  with  increasing  plastic  deformation,  there  is  a  considerable  variation  of  J  through  the  thick¬ 
ness  with  the  value  at  the  center  being  much  higher  than  that  at  the  free  surface.  Similar  trends 
have  been  reported  in  Nakamura,  Shih  and  Freund  (1986a)  and  Delorenzi  and  Shih  (1983). 
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43  Crack  Tip  Opening  Displacement 

The  crack  tip  opening  displacement  was  calculated  from  the  numerical  solution  using  the 
45°  intercept  procedure  introduced  by  Shih  (1981).  The  crack  tip  opening  displacement  at  the 
centerline  of  the  3-D  specimen  is  plotted  against  the  load  in  Fig.  5.  The  plane  stress  and  plane 
strain  results  are  also  shown  in  the  figure  for  comparison.  As  before,  the  3-D  curve  is  bounded 
by  the  2-D  analyses. 

The  variation  of  the  crack  tip  opening  displacement  along  the  crack  front  is  shown  in  Fig.  6 
for  three  different  load  ratios  P/Pe.  The  crack  tip  opening  displacement,  5,  has  been  made 
dimensionless  by  the  value  5Ci,  at  the  centerline  of  the  specimen.  As  in  the  case  of  J  in  Fig. 
4,  There  is  a  strong  variation  of  the  crack  tip  opening  with  increasing  plastic  deformation,  with 
the  value  at  the  centerline  being  very  much  in  excess  of  that  at  the  free  surface.  Such  strong 
variation  of  crack  tip  opening  through  the  thickness  was  also  observed  in  Delorenzi  and  Shih 
(1983)  and  Wellman,  Rolfe  and  Dodds  (1987).  This  result  along  with  Fig.  4  suggests  the  possi¬ 
bility  of  a  tunneling  mode  of  fracture  with  the  crack  beginning  to  propagate  first  at  the  center. 
Indeed,  some  evidence  of  this  could  be  seen  by  a  post  mortem  examination  of  the  fracture  sur¬ 
face  of  the  experimental  specimen.  In  Fig.  7,  the  crack  opening  displacement  at  the  center  of 
the  specimen,  which  is  normalized  by  (J/a0)  is  plotted  against  the  quantity  J/(a0C).  Here  C  is 
the  uncracked  ligament  length.  The  2-D  results  are  also  shown  in  the  figure.  The  dimension¬ 
less  quantity  J/(o  „  C )  is  an  indication  of  the  extend  of  plastic  yielding  (Shih,  1981  and  Shih  and 
German,  1981).  Small  values  of  J/(ocC )  (less  than  about  0.001)  indicate  well-contained  yield¬ 
ing,  while  large  values  (greater  than  0.01)  correspond  to  large  scale  plasticity. 

The  values  of  5,/(J /o0)  based  on  the  2-D  solutions  are  fairly  constant  over  the  entire  range 
of  plastic  deformation.  For  plane  stress,  8,/(J/c0)  varies  from  0.82  for  well  contained  yielding 
to  about  0.80  for  large  scale  plasticity,  while  for  plane  strain  it  varies  from  0.62  to  about  0.60. 
These  values  are  somewhat  smaller  than  the  HRR  solution  which  gives  0.95  and  0.68  for  plane 
stress  and  plane  strain  respectively  corresponding  to  n  =  22  (Shih,  1981  and  Shih,  1983). 

However,  Narasimhan  and  Rosakis  (1988),  reported  a  value  of  8,/(J/o0)  of  0.85  for  the 
non-hardening  case  based  on  a  finite  element  solution  which  modelled  plane  stress,  small  scale 
yielding  conditions.  Shih  (1981)  has  also  computed  the  crack  opening  displacement  for  dif¬ 
ferent  hardening  exponents  and  specimen  geometries  based  on  plane  strain  finite  element 
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solutions.  For  the  Cracked  Bend  Bar  (of  perfectly  plastic  material),  he  reports  a  value  of 
6,/(J/o0)  which  varies  from  0.65  to  0.58  from  contained  yielding  to  fully  plastic  conditions. 
As  noted  in  Narasimhan  and  Rosakis  (1988),  the  discrepancy  with  the  HRR  solution  in  the 
non-hardening  limit  is  due  to  the  fact  that  the  HRR  strain  solution  may  not  dominate  over  any 
finite  distance  near  the  tip. 

It  can  be  seen  from  Fig.  7  that  the  normalized  crack  opening  displacement,  5,/(J/a0)  at  the 
centerline  of  the  3-D  specimen  varies  from  about  0.68  for  contained  yielding  to  about  0.75  for 
large  scale  yielding.  For  contained  yielding,  this  value  is  closer  to  plane  strain  whereas  for 
large  scale  yielding  is  closer  to  the  plane  stress  limit.  A  qualitatively  similar  trend  may  be 
noticed  in  the  results  obtained  in  Delorenzi  and  Shih  (1983)  for  a  compact  tension  specimen. 

4.4  Stress  Variation  Through  the  Thickness 

The  variation  of  the  plane  strain  constraint  (as  given  by  the  ratio  o33/(gu  +  o  22)  through 
the  thickness  of  the  specimen  is  shown  in  Fig.  8  for  different  load  ratios  P  /  PQ.  The  distribu¬ 
tion  for  the  elastic  case  is  also  plotted  in  the  figure.  The  stresses  in  this  plot  were  evaluated  at 
the  centroid  of  the  elements  nearest  to  the  tip,  located  at  r/h  =  0.005  and  0  =  0°. 

The  first  observation  that  can  be  made  is  that  the  ratio  a  33/(0  n  +  022)  is  uniform  over  most 
of  the  specimen  thickness  (along  the  crack  front)  and  drops  to  zero  very  near  the  free  surface. 
Its  value  in  the  interior  of  the  specimen  is  about  0.28  for  the  elastic  case  and  approaches  0.47 
with  increasing  plastic  deformation.  It  should  be  noted  here  that  under  plane  strain  conditions 
o  33/(0  „  +  O22)  is  equal  to  the  Poissons  ratio  (which  is  0.3  in  the  present  analysis)  for  the  elas¬ 
tic  case.  Also  this  ratio  is  equal  to  0.5  in  the  fully  plastic  limit. 

Secondly,  the  numerical  results  indicate  that  even  for  an  applied  load  of  only  0.45  P0  the 
ratio  0  33/(  On  +  O22)  near  the  crack  front  is  0.46  in  the  interior  of  the  specimen.  The  ratio 
shows  very  little  variation  with  increase  in  load  level.  This  is  due  to  the  low  hardening  of  the 
material  which  allows  very  large  in-plane  plastic  strains  to  accumulate  near  the  tip  even  for 
moderate  extents  of  yielding.  This  leads  to  a  stress  distribution  near  the  crack  front  in  the  inte¬ 
rior  of  the  specimen  which  is  similar  to  the  plane  strain  Prandtl  field,  (see  Rice  and  Rosengren, 
1968,  and  Rice,  1968),  also  see  Sec.  4  and  Fig  22.  The  observation  that  c33  /  (  On  +  o 22) 
approaches  the  value  0.5  near  the  crack  front,  is  consistent  with  the  expectation  of  a  plane  strain 
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like  constraint  in  the  center  plane  of  the  specimen. 

The  variation  through  the  thickness  of  the  cartesian  stress  components  On.  On  and  o33, 
normalized  by  the  yield  stress  a is  shown  in  Figs.  9,  10,  and  1 1  for  three  different  load  levels 
ranging  from  contained  yielding  to  large  scale  plasticity.  The  stresses  displayed  in  these 
figures  are  computed  at  0  =  0  0  and  r/h  =  0.005,  0.025,  0.065,  0.125,  0.205,  0.305,  0.425  and 
0.565.  The  arrow  direction  in  the  figures  indicates  increasing  r/h  values.  The  distribution 
through  the  thickness  of  the  normalized  triaxial  stress  (o  n  +  o22  +  a33)  /  3  oa  is  also  shown 
in  these  figures. 

It  can  be  observed  that  very  near  the  crack  front  (in  the  interior  of  the  specimen)  the  stresses 
are  highly  elevated  above  the  yield  stress  leading  to  a  strongly  triaxial  stress  field.  This  is 
characteristic  of  the  plane  strain  crack  tip  fields  (Rice  and  Rosengren,  1968,  and  Rice,  1968). 
For  low  load  levels  (Fig.  9),  the  stress  variation  along  the  crack  front  is  uniform  over  most  of 
the  specimen  thickness  and  drops  near  the  free  surface.  But  in  the  fully  plastic  limit  (Fig.  1 1)  a 
strong  variation  of  stresses  along  the  crack  front  is  noticed.  This  leads  to  a  considerable  loss  of 
triaxiality  as  the  free  surface  is  approached  along  the  crack  front.  This  observation  is  important 
since  a  strongly  triaxial  field  influences  the  growth  and  nucleation  of  voids. 

These  figures  show  that  for  r/h  >  0.5,  a  33  is  close  to  zero  uniformly  over  the  specimen 
thickness  irrespective  of  the  load  level.  Also,  the  other  stress  components  have  very  little  varia¬ 
tion  through  the  thickness  for  r/h  >  0.5.  This  suggests  that  plane  stress  conditions  are 
approached  for  distances  from  the  crack  front  in  excess  of  half  the  specimen  thickness  (see  Sec. 
5  for  further  discussion). 

4.5  Surface  (out-of-plane)  displacements 

Two  views  of  the  surface  displacement  field  u3  of  the  deformed  specimen  at  a  load  of 
52300  N  (obtained  from  the  numerical  calculation)  are  shown  in  Fig.  12.  The  out-of-plane  sur¬ 
face  displacements  of  Fig.  12  were  used  to  generate  numerical  interferograms  for  comparison 
with  the  experimental  results  described  in  Part  II.  Figure  13  is  a  composite  picture  showing 
one-half  of  an  experimental  and  one-half  of  a  numerically  generated  interferogram.  As  is  evi¬ 
dent  from  this  comparison,  good  qualitative  agreement  of  experimental  and  synthetic  interfero¬ 
grams  is  achieved.  In  Fig.  14,  the  numerical  and  expenmental  out-of-plane  displacements  along 
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the  ray  ahead  of  the  cracktip,  u  3  (  r,  9  =  0),  are  plotted  against  normalized  distance  from  the  tip 
at  a  load  of  52300  N.  The  displacements  have  been  made  dimensionless  by  the  quantity 
(J/o0).  The  extent  of  the  plastic  zone  ahead  of  the  tip  as  well  as  the  location  along  the  x j  axis 
corresponding  to  half  the  specimen  thickness  are  indicated  for  reference.  The  out-of¬ 
displacement  corresponding  to  the  plane  stress  HRR  field  at  the  same  average  J  value  is  also 
plotted  for  comparison. 

It  can  be  observed  that  there  is  excellent  agreement  between  the  3-D  numerical  and  experi¬ 
mental  results,  even  very  near  the  crack  tip.  The  agreement  was  found  to  be  very  good  over  the 
entire  range  of  loading  (see  part  II).  It  was  found  that  the  plane  stress  HRR  field  did  not  ade¬ 
quately  describe  u  3  on  the  specimen  surface.  The  above  observation  has  also  been  made  dur¬ 
ing  previous  experimental  investigations  using  the  optical  method  of  caustics  (Rosakis,  Zehnder 
and  Narasimhan,  1987  and  Zehnder,  1987).  In  these  studies,  it  was  noticed  that  caustics  based 
on  plane  stress  HRR  fields  could  not  always  be  used  to  accurately  measure  /.  This  discrepancy 
is  because  of  three  dimensionality  near  the  crack  tip  and  due  to  finite  specimen  dimensions  that 
cause  higher  order  terms  to  become  important  away  from  the  crack  tip.  Further  comparisons 
between  numerical  and  experimental  results  based  on  caustics  and  interferometry  may  be  found 
in  Part  II. 


5.  THREE  DIMENSIONAL  EFFECTS  AND  HRR  DOMINANCE 

5.1  Radial  Stress  Variation 

The  radial  variation  of  the  normalized  opening  stress,  cr  22  /  e  >  along  the  uncracked  liga¬ 
ment  obtained  from  the  plane  stress  and  plane  strain  numerical  solutions,  are  displayed  in  Figs. 
15  and  16  respectively.  Results  are  given  for  three  load  levels  ranging  from  well  contained 
yielding  to  large  scale  plasticity.  The  extent  of  plastic  deformation  is  indicated  by  the  parame¬ 
ter  C<j0/J  (see  Sec.  4.3).  The  stress  distribution  based  on  the  HRR  field  is  also  shown  in  these 
figures. 

It  can  be  seen  from  Fig.  15  that  in  the  case  of  plane  stress,  the  HRR  stress  field  agrees  very 
well  with  the  numerical  solution  over  a  large  distance  (x  1  less  than  about  30  J/80)  for  all  levels 
of  plastic  deformation.  A  similar  observation  was  made  in  Narasimhan  and  Rosakis  (1988)  for 
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plane  stress,  small  scale  yielding  conditions.  Hence  it  appears  that  under  plane  stress,  HRR 
dominance  is  assured  for  distances  along  the  ligament  of  the  order  of  30  to  40  times  the  crack 
tip  opening  displacement  6, . 

The  radial  stress  distribution  for  plane  strain  (Fig.  16)  suggests  that  the  HRR  field  dominates 
only  very  near  the  tip  (for  x  \  less  than  about  3  J/cD).  This  is  in  good  agreement  with  the 
results  of  Shih  and  German  (1981)  for  the  Cracked  Bend  Bar.  Thus,  as  noted  in  Shih  and  Ger¬ 
man  (1981),  HRR  dominance  to  about  6  to  10  8,  for  all  levels  of  plastic  yielding. 

The  radial  variation  of  the  three  dimensional  normalized  out-of-plane  stress  component 
<*33  /  (  <*  it  +  O22)  (which  is  an  indication  of  the  plane  strain  constraint)  ahead  of  the  crack  front 
is  shown  in  Fig.  17  along  three  different  planes  through  the  specimen  thickness.  These  planes 
correspond  to  x3  /h  of  0.075  (which  is  near  the  mid  plane  of  the  specimen),  0.325  (near  the 
quarter  plane)  and  0.475  (near  the  free  surface).  These  results  have  been  taken  from  the  3-D 
numerical  solution  at  applied  load  levels  of  0.45P„,  0.7Po  and  P0. 

As  noted  earlier(  in  connection  with  Fig.  8),  it  can  be  seen  from  this  figure  that  the  ratio 
a  33/(0  ii  +  a  22)  approaches  a  value  close  to  0.5  near  the  crack  tip  over  most  of  the  specimen 
thickness  once  moderate  amounts  of  yielding  take  place.  Also,  irrespective  of  load  level  and  for 
all  planes  through  the  thickness,  the  above  ratio  becomes  quite  small  for  radial  distances  from 
the  crack  tip  exceeding  approximately  half  the  specimen  thickness  (r  >  0.5  h).  This  observation 
along  with  Figs.  9,  10  and  11,  which  show  that  the  other  stress  components  have  negligible 
variation  through  the  thickness  for  r/h  >  0.5,  suggests  that  plane  stress  conditions  are 
approached  for  distances  from  the  crack  tip  exceeding  half  the  specimen  thickness.  A  similar 
conclusion  was  reached  by  Yang  and  Freund  (1984)  on  the  basis  of  an  analytical  investigation 
and  by  Parsons,  Hall  and  Rosakis  (1986),  on  the  basis  of  a  numerical  study  of  3-D,  elastic 
cracked  specimens.  Experimental  observations  based  on  the  optical  method  of  caustics 
(Rosakis  and  Ravichandar,  1984,  and  Rosakis,  Zehnder  and  Narasimhan,  1987),  also  indicate 
that  three  dimensional  effects  become  important  only  for  r/h  <  0.5. 

In  order  to  examine  the  above  issue  further,  the  radial  distribution  of  the  normalized  opening 
stress  O22  /  <*<7  ahead  of  the  crack  front  is  shown  along  three  different  planes  through  the  speci¬ 
men  thickness  x3  /  h  =  0.075,  0.325  and  0.475)  in  Fig.  18.  The  bottom  scale  in  the  figure  gives 
the  radial  distance  from  the  tip  normalized  by  the  specimen  thickness  h,  while  in  the  top  scale 
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the  nonnalization  is  done  with  respect  to  the  uncracked  ligament  length  C.  For  comparison,  the 
corresponding  variation  based  on  the  full  field  plane  stress  and  the  full  field  plane  strain  numeri¬ 
cal  solutions  (at  the  same  load  level  P)  are  indicated  in  the  figure.  Results  are  once  again  pro¬ 
vided  for  three  values  of  applied  load. 

This  figure  shows  that  for  r/h  >  0.5,  the  plane  stress  radial  stress  distribution  is  approached 
over  the  entire  specimen  thickness  irrespective  of  the  load.  On  the  other  hand,  the  plane  strain 
radial  stress  variation  provides  a  good  approximation  near  the  mid-plane  of  the  specimen  for  r/h 
<  0.2.  It  should  be  noted  that  the  o  22  (opening  stress)  distribution  is  highly  elevated  above  the 
yield  stress  as  the  crack  tip  is  approached  in  the  interior  of  the  specimen.  Also,  it  can  be  seen 
that  there  is  substantial  variation  of  stresses  through  the  thickness  as  fully  plastic  conditions  are 
approached  (see  Fig.  11  and  distribution  corresponding  to  P  =  P„  in  Fig.  18). 

In  Fig.  19,  the  radial  variation  of  the  opening  stress  ahead  of  the  crack  tip  is  again  shown  in 
the  same  non-dimensional  scales  but  is  now  compared  with  the  analytical  asymptotic  plane 
stress  and  plane  strain  HRR  distributions.  The  asymptotic  (HRR)  distributions  corresponding  to 
the  same  average  value  of  J  have  been  plotted  in  this  figure.  This  figure  demonstrates  that  very 
near  the  crack  tip  (r  <  0.02  h)  the  opening  stress  in  the  interior  of  the  specimen  tends  towards 
the  asymptotic  plane  strain  HRR  solution.  On  the  other  hand,  the  plane  stress  HRR  solution  is 
approached  away  from  the  crack  tip  within  the  plastic  zone.  This  is  more  evident  with  increas¬ 
ing  plastic  yielding  (see  figures  corresponding  to  P  =  0.7Po  and  P  =  Pa). 

S2  Angular  stress  variation 

The  angular  variation  of  stresses  in  the  ring  of  elements  closest  to  the  crack  tip 
(r  =0.005 h  =0.001  C)  obtained  from  the  plane  stress  and  plane  strain  numerical  solutions  is 
displayed  in  Figs.  21  and  22  respectively.  The  stresses  in  these  plots  have  been  normalized  by 

-  3  >A 

the  maximum  value  of  the  equivalent  stress  a  =  (ySySy)  .  Th  angular  stress  variation  given 

by  the  HRR  solution  is  also  shown  in  these  figures  for  comparison.  The  values  of  the  load 
ratios  P IP 0  and  the  dimensionless  quantity  r/(J !o0)  are  also  indicated  in  the  figures.  It  can  be 
seen  that  the  near-tip  two-dimensional  numerical  results  display  very  good  agreement  with  the 
corresponding  HRR  distribution  for  the  entire  loading  range. 
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In  Fig.  22,  the  near-tip  angular  distribution  of  the  normalized  stresses  obtained  from  the  3-D 
numerical  solution  at  *3  =  0.075/:  (near  the  mid-plane  of  the  specimen)  and  r  =  0.005  h  is 
shown.  Results  are  given  for  three  different  load  ratios  P IP0 .  Also  plotted  in  the  figure  is  the 
angular  stress  variation  given  by  the  plane  strain  HRR  asymptotic  solution.  It  is  found  that  the 
3-D  crack  tip  stresses  near  the  mid-plane  of  the  specimen  agree  very  closely  with  the  plane 
strain  HRR  field. 

It  should  be  recalled  that  the  material  used  in  the  experiments,  and  modeled  by  the  numeri¬ 
cal  calculations,  has  a  very  low  strain  hardening  (n  =22)  and  hance  the  angular  stress  variation 
in  Fig.  22  displays  all  the  essential  features  of  the  perfectly-plastic  plane-strain  Prandtl  slip  line 
field.  For  example,  in  the  angular  range  45°  <  0  <  135°,  <sn  =  Oee  and  or6  is  constant.  These 
are  characteristic  of  the  "centered  fan"  region  of  the  Prandtl  field  (Rice,  1968)  where  the  plastic 
strains  are  singular  as  1/r.  Also,  note  that  the  stresses  are  highly  elevated  (above  amiX)  ahead 
of  the  crack  tip  (9  =  0°).  This  result  clearly  demonstrates  the  dominance  of  the  plane  strain 
HRR  field  along  the  crack  front  in  the  interior  of  the  specimen. 

S3  Plastic  zones 

The  plastic  zones  corresponding  to  the  two  different  load  levels  and  three  different  planes 
through  the  specimen  thickness  (xj/h  =  0.075,0.325  and  0.475)  are  shown  in  Figures  23  and  24. 

Both  figures  display  contours  of  equal  equivalent  stress  a.  The  contour  further  away  from 
the  crack  tip  (boundary  of  the  zones)  corresponds  to  0  =  0.99  o„ .  The  remaining  contours 
correspond  to  <j=1.02cto,  1.05  o0,  1.08  oc,  1.11  o0 ,  respectively.  Figure  23  corresponds  to 
P  -  0.7  Pa .  For  this  load  level,  the  plastic  zone  has  not  yet  interracted  with  the  applied  loads  at 
the  boundary.  For  all  cross  sections,  the  maximum  plastic  zone  site  is  less  than  half  the  speci¬ 
men  thickness  and  thus  it  lies  well  within  the  region  of  high  near  tip  three  dimensionality  as 
discussed  in  section  5.1.  It  is  somewhat  surprising  to  note  here  that  the  plastic  zone  shape 
obtained  at  the  specimen  surface  looks  closer  to  the  plane  strain  plastic  zone  shape  than  its 
counterpart  in  the  middle.  This  contradicts  the  popular  concept  of  plastic  zone  shapes  as 
reported  in  a  number  of  books  on  fracture  mechanics,  (Hellan,  1985,  Kanninen  and  Popelar, 
1985). 
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Figure  24  shows  contours  of  o  for  P  -  Pa.  The  figure  clearly  indicates  the  interractions  of 
the  crack  tip  plastic  zone  with  the  plastic  deformation  at  the  specimen  boundary  due  to  the 
applied  loads.  In  this  case  and  for  all  cross  sections,  the  outer  contour  of  equal  c  reveals  almost 
circular  plastic  zone  shapes  reminiscent  of  the  shapes  predicted  by  Narasimhan  and  Rosakis 
(1988),  on  the  basis  of  a  plane  stress  analysis.  This  is  not  surprising  since  the  elastic  plastic 
boundary  now  lies  well  outside  a  circle  of  radius  hll  (dotted  line)  which  signifies  the  region 
outside  which  plane  stress  conditions  dominate.  Indeed,  as  shown  in  figure  17,  the  stress  field 
for  all  load  levels,  approaches  the  plane  stress  constraint  for  rlh  >  0.5. 

6.  CONCLUSIONS 

A  simultaneous  experimental  and  numerical  (both  3-D  and  2-D)  investigation  has  been  per¬ 
formed  to  simulate  the  monotonic  loading  of  a  ductile  3-point  bend  specimen  from  contained  to 
fully  plastic  conditions.  The  average  J  and  the  crack  tip  opening  displacement  from  the  3-D 
analysis  are  bounded  by  the  2-D  (plane  stress  and  olane  strain)  analyses.  The  cracktip  opening 
displacement  and  stresses  exhibit  a  strong  variation  along  the  crack  front  with  increasing  plastic 
deformation.  The  3-D  numerical  results  clearly  indicate  that  the  plane  strain  HRR  field  dom¬ 
inates  very  near  the  crack  front  (r  <  0.02  h )  in  the  interior  of  the  specimen.  On  the  other  hand, 
the  plane  stress  HRR  field  agrees  well  with  the  3-D  numerical  results  obtained  on  the  surface 
layer  of  elements.  The  near-tip  angular  stress  variation  close  to  the  mid-surface  of  the  specimen 
displays  all  the  essential  features  of  the  asymptotic  plane  strain  Prandtl  field,  once  moderate 
amount  of  yielding  has  occurred. 

On  the  other  hand,  there  is  a  considerable  loss  in  triaxiality  of  the  stress  field  as  the  free  sur¬ 
face  is  apporached  along  the  crack  front.  This  is  parliculary  evident  for  the  larger  applied  loads 
and  suggests  a  tendancy  towards  loss  of  local  J  dominance  along  the  crack  front  as  large  scale 
yielding  conditions  are  approached.  The  radial  and  thickness  variations  of  the  stresses  from  the 
3-D  analysis  suggests  that  for  r  >0.5  h,  plane  stress  conditions  are  approached.  Comparison  of 
the  results  to  2-D  analyses,  whenever  appropriate,  and  the  excellent  agreement  between  the 
numerical  and  experimental,  P  -6,  J  -P  and  Mj(r ,  8)  results  leave  no  question  that  the  3-D 
calculations  are  accurate.  These  comparisons  also  show  that  the  calculation  is  sufficiently 
refined  to  model  this  specimen  well,  both  in  the  overall  behavior  and  in  the  details  of  the  near¬ 
tip  fields. 
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Figure  2.:  Load  -  Load  displacement  curve. 


Plane  Stress 


Figure  3a.:  J  versus  applied  load 


versus  load  point  displacement. 


Figure  4.:  Variation  of  j  through  the  thickness  for  different  load  ratios.  P0  is  the  plane 

stress  limit  load  for  a  rigid-perfectly  plastic  material. 
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Figure  5.:  Crack  tip  opening  displacement  at  centerline  versus  load. 


Figure  6.:  Variation  of  the  crack  opening  displacement  along  the  crack  front. 


normalized  centerline  crack  tip  opening  displacement 
rielding  to  fully  plastic  conditions. 


P/P.  •  0. 7 


Variation  of  die  plane  strain  constraint,  o33  /  (a 


Variation  of  normalized  stress  components  through  the  thickness  for  load  ratio 
P/P0  =  0.45.  The  curves  shown  correspond  to  r/h  -  0.005,  0.025,  0.065,  0.125, 
0.205,  0.305,  0.425  and  0.565.  The  arrow  direction  indicates  increasing  r/h 


P/Pc*  0.7 


Variation  of  normalized  stress  components  through  the  thickness  for  load  ratio 
P/P o  =  0.7.  The  curves  shown  correspond  to  r/h  -  0.005,  0.025,  0.065,  0.125, 
0.205,  0.305,  0.425  and  0.565.  The  arrow  direction  indicates  increasing  r/h 


d/d 


Variation  of  normalized  stress  components  through  the  thickness  for  load  ratio 
P/P0  =  1.0.  The  curves  shown  correspond  to  r/h  -  0.005,  0.025,  0.065,  0.125, 
0.205,  0.305,  0.425  and  0.565.  The  arrow  direction  indicates  increasing  r/h 


Figure  12a.:  Out-of-plane  displacement  u  3  for  load  of  52300  N. 


32300 


Figure  14.:  Nondimensional  u3  displacement  on  the  0  =  0  line.  Comparison  of  experimental 
and  3-D  numerical  results  for  load  of  52300  N.  Also  shown  is  the  HRR  field 
within  the  plastic  zone. 
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versus  normalized  distance  ahead  of  the  crack  tip.  Comparison  between 
rain  numerical  and  HRR  solutions  from  contained  to  fully  plastic  condi- 
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variation  of  plane  strain  constraint  a33/(an  +  o22)  along  three  different 
through  the  specimen  thickness.  Results  are  given  for  three  different  load 


ion  of  o  22  /  a  „  along  three  different 


P/p0  =  0.45 


Radial  variation  of  a  22  /  oQ  along  three  different  planes  through  the  specimen 
thickness.  Comparison  between  3-D  numerical  and  2-D  HRR  solutions. 


angular  stress  variation.  Comparison  between  plane  stress  numerical 
solutions  from  contained  to  fully  plastic  conditions. 
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Comparison  between  3-D  near-tip  angular  stress  distribution  along  the  plane 
x  3  /  h  =  0.075  and  HRR  plane  strain  solution. 
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